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Abstract 

Electromagnetic wave scattering by many parallel to z—axis, thin, 
perfectly conducting, circular infinite cylinders is studied asymptotically 
as a — > 0. Let Dm be the crossection of the m— th cylinder, a be its 
radius, and Xm — i^Xmi^ Xm,2^ b6 its C6nt6r, 1 ^ 771 ^ ^ — A^(ci^. It is 
assumed that the points Xm are distributed so that 

A/'(A) = lni [ N{x)dx[l + o{l)], 

where A/'(A) is the number of points in an arbitrary open subset A 
of the plane xoy. The function N{x) > is a given continuous function. 
An equation for the self-consistent (efficient) field is derived as a 0. 
The cylinders are assumed perfectly conducting. A formula is derived for 
the effective refraction coefficient in the medium in which many cylinders 
are distributed. These cylinders may model nanowires embedded in the 
medium. Our result shows how these cylinders influence the refraction 
coefficient of the medium. 

PACS: 03.40.Kf, 03.50 De, 41.20. Jb, 71M.+C 
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1 Introduction 

There is a large fiterature on electromagnetic (EM) wave scattering by an array 
of parallel cylinders (see, e.g., [2], where there are many references given, and 
[3]). Electromagnetic wave scattering by many parallel to z—axis, thin, per- 
fectly conducting, circular, of radius a, infinite cylinders is studied in this paper 
asymptotically as a -> 0. The cylinders are thin in the sense ka ^ 1, where k 
is the wave number in the exterior of the cylinders, 
The novel points in this paper include: 

1) The solution to the wave scattering problem is considered in the limit 
a —7- when the number M = M{a) of the cylinders tends to infinity at a 
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suitable rate. The equation for the limiting (as a — ?> 0) effective (self-consistent) 
field in the medium is derived, 

2) This theory is a basis for a method for changing refraction coefficient in 
a medium. The thin cylinders model nanowires embedded in the medium. The 
basic physical resul of this paper is formula (j48p , which shows how the embedded 
thin cylinders change the refraction coefficient n'^{x). 

Some extension of the author's results ([7]- [IS]) is obtained for EM wave 
scattering by many thin perfectly conducting cylinders. 

Let Dm, 1 < m < M, be a set of non- intersecting domains on a plane P, 
which is xoy plane. Let Xm G Dm, Xm — {xmi, Xm2), be a point inside and 
Cm be the cylinder with the cross-section Dm and the axis, parallel to z-axis, 
passing through Xm ■ We assume that Xm is the center of the disc Dm if Dm is 
a disc of radius a. 

Let us assume that the cylinders are perfect conductors. Let a = O.SdiamZ?™. 
Our basic assumptions are 

A:a<l, (1) 

where k is the wave number in the region exterior to the union of the cylinders, 
and 

AA(A)=ln-/ N{x)dx[l + o{l)], a ^ 0, (2) 
Ja 

where A/'(A) = 1 is the number of the cylinders in an arbitrary open 

subset of the plane P, N{x) > is a continuous function, which can be chosen as 
we wish. The points distributed in an arbitrary large but fixed bounded 

domain on the plane P. We denote by the union of domains Dm, by ft' its 
complement in P, and by D' the complement of D in P. The complement in 
of the union C of the cylinders Cm we denote by C. 

The EM wave scattering problem consists of finding the solution to Maxwell's 
equations 

V X E = iuj^iH, (3) 
y X H = -iuoeE, (4) 

in C , such that 

Et^O on dC, (5) 

where dC is the union of the surfaces of the cylinders Cm , Et is the tangential 
component of E, fj, and e are constants in C", w is the frequency, fc^ — uj^^^jl, k 
is the wave number. Denote by tiq = e/i, so fc^ — w^rtg. The solution to ©-(O 
must have the following form 

E{x) ^ E(i[x) +v(x), X ^ {xi,X2,X3) ^ {x,y,z) ^ {x,z), (6) 

where Eq{x) is the incident field, and v is the scattered field satisfying the 
radiation condition 

V~r{^^^-^kv^^o{l), r^{x\^x\fl\ (7) 
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and we assume that 



{cj}, j ~ 1, 2, 3, are the unit vectors along the Cartesian coordinate axes a:, y, z. 
We consider EM waves with := = 0, i.e., E-waves, or TH waves, 

-A V X 

E = YEjej, H = Hiei + H2e2 ^ . (9) 

One can prove (see Appendix) that the components of E can be expressed by 
the formulas: 

Ej^'-^U^^e^'^'^-, J = 1,2, £;3 = C/e^'=^^ U = ju, (10) 
where '■= ■§^, u — u{x, y) solves the problem 

(As + n^)u = in 17' (11) 

"Ian = 0, (12) 

u = e'^y + (13) 

and w satisfies the radiation condition ([7]). Similar calculations are done with 
Borgnis potentials (see, e.g., |1]). The unique solution to ©-(jH]) is given by the 
formulas: 

i?i = ^C/.e^'^^ E^^'-^Uye''^'^ i?3 = C/e*'=^^ (14) 

i?i = -^C/,e'^■^^ = — ^C/xe*'■^^ i/3=0, (15) 

where Ux '■— Uy is defined similarly, and u = u{x) = u{x,y) solves scalar 
two-dimensional problem ([TT]) -([T ^ . These formulas are derived in the Appendix 
for convenience of the reader. 

Problem (fTT|) - p3)) has a unique solution (see, e.g., [4]). 

Our goal is to derive an asymptotic formula for this solution as a — > 0. Our 
results include formulas for the solution to the scattering problem, derivation of 
the equation for the effective field in the medium obtained by embedding many 
thin perfectly conducting cylinders, and a formula for the refraction coefficient in 
this limiting medium. This formula shows that by choosing suitable distribution 
of the cylinders, one can change the refraction coefficient, one can make it 
smaller than the original one. 

The paper is organized as follows. 

In Section 2 we derive an asymptotic formula for the solution to (ITT]) - P^ 
when M = 1, i.e., for scattering by one cylinder. 

In Section 3 we derive a linear algebraic system for finding some numbers 
that define the solution to problem ([TT1) -(fT3 |) with M > 1. Also in Section 3 we 
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derive an integral equation for the effective (self-consistent) field in the medium 
with M{a) — >■ oo cyhnders as a — >■ 0. At the end of Section 3 these results are 
applied to the problem of changing the refraction coefficient of a given material 
by embedding many thin perfectly conducting cylinders into it. 

In Section 4 conclusions are formulated. 

In Appendix formulas p^ - p^ are derived. 

2 EM wave scattering by one thin perfectly con- 
ducting cylinder 

Consider problem (ITT]) - p^ with ft = Di, il' being the complement to Di in 
R^. Assume for simplicity that Di is a circle xf + X2 < a^- 
Let us look for a solution of the form 

g{x,t)a{t)dt, g{x,t):=^H^'\K\x-t\), (16) 
where Si is the boundary of Di, H^^^ is the Hankel function of order 1, with 



index 0, and a is to be found from the boundary condition ()12p . It is known 
(see, e.g., [S]) that 

^(Kr) = a(K) + -!-ln- + 0(1), as r ^ 0, (17) 
2n r 

where 

a(«):=^ + ;lln-. (18) 

The function satisfies equations pT|) and for any a, and if a is such that 
function pB|) satisfies boundary condition p^ . then u solves problem pT|) - p^ . 
We assume a sufficiently smooth (Holder-continuous is sufficient). 

The solution to problem pTjl -pB ]) is known to be unique (see, e.g., [4]). 
Boundary condition p2|) yields 



- uois) = aiK)Q + go{s,t)a{t)dt, Q := / cT(t)dt, (19) 
Jsi Jsi 

uo(s) :=e*""^ se go{s,t) = ^\n—, rst '-^ \s - t\. (20) 

27r rst 

If ka < 1, /c^ = + fc|, then 

uo{s) = 1 + 0{Ka). 

Equation p9)) is uniquely solvable for tr if a is sufficiently small [5]. 
We are interested in finding asymptotics of Q as a — >■ 0, because u{x) in pB|) 
can be well approximated in the region |a;| 3> a by the formula 

u{x) ^ U(^{x) + g{x,Q)Q + o{l), a 0, 
Uo{x)^e"'''\ X2 = y. 
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To find asymptotics of Q as a 0, let us integrate equation over 5*1 and 
obtain 



f dt(j{t)^ [ \nrstds, (22) 

J Si J Si 



where IS*!! is the length of 5*1, l^il = 27ra if 5*1 is the circle |a;| = a, and 
rst = \s — t\. Denote 



/ := — / Inrstdt = 0{alna), a^O. 
Jsi 

If Si is the circle |a;| = a, integral ([^ can be calculated analytically: 

/ = — / In ^2a^ ~ 2a^ cos(V' — (p)dip 
27r Jo 



(23) 



2-!T r2TT I I 



— / \ny2a?dif+— j In y 2 sin — - — dip. 
Ztt ,Iq ztt 



(24) 



Thus, 



/-aln(V2a) + ^ + — In | sin (25) 
2 27r ,/n 2 



One can derive that 

/ ln|sin— = 2 / In|sin6'|d6' = -27rln2. (26) 
Indeed, if J := /J^ In | sin9\d0, then 



J= / ln2|sin-cos-|(i6' = 7rln2+ / lnsin-d6l+/ In | cos -|d6l = 7rln2+2J. 
Jo 2 2 7q 2 Jo 2 

Thus, J = -7rln2. From ^ and ^ one gets 

/ = alna(l + o(-^ ) , a ^ 0. (27) 



From (111) and (HH) it follows that 



In a\ 



Therefore, the asymptotic solution to the scattering problem (|lip - p^ in the 
case of one circular cylinder of radius a, as a — >■ 0, is 

, ^ 27r , , , , 
u{x) uo{x) - —^-g{x,0)uo{0), a 0, > a. (29) 
in - 

a 

Electromagnetic wave, scattered by the single cylinder, is calculated by formulas 
(fT4 |) -(fT5 |) in which u — u{x) u{xi, X2) is given by formula (|29p . 
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3 Wave scattering by many thin cylinders 

Problem (Ilip - (I13I) should be solved when $7 is a union of many small domains 
Dm, ^ = ^m=iDm- We assume that Z?,„ is a circle of radius a centered at the 
point Xm- 

Let us look for u of the form 

M „ 

u{x) ^ un{x) + / g{x,t)a„i{t)dt. (30) 

m=l •^■Sm 

We assume that the points distributed in a bounded domain D on the 

plane P = xoy by formula ([2]). The field uq{x) is the same as in Section 2, 
uo(x) = e*''^, and Green's function g is the same as in formulas ([TC|) -([T ^ . It 
follows from © that M = M{a) = 0{\n\). We define the effective field, acting 
on the Dj by the formula 



Si 



M 



Ue — u'^J^ — u{x) — I g{x,t)cjj{t)dt, \x — Xj \ > a, (31) 

which can also be written as 

Ue{x) = UQ{x)+ y2 / 9ii,t)crmit)dt. 

/ • J s„. 

We assume that the distance d = d{a) between neighboring cylinders is much 
greater than a: 

d:S>a, lim = 0. (32) 

Let us rewrite ([SO]) as 

M M 



U — Uq + 



g{x,Xm)Q„i+ / [g{x,t) - g{x,Xm)]cr7n{t)dt, (33) 

, 1 ™ I-^Stti 



m— 1 m— 1 



where 

Qrn / <Jm{t)dt. (34) 

As a 0, the second sum in (|33)) (let us denote it i;2) is negligible compared 
with the first sum in (j33p . denoted Ei, 

iSal^lSil, a^O. (35) 

The proof of this is similar to the one given in [5] for a similar problem in R'^. 
Let us check that 

\gix,Xm)Qm\^ \ [gix^t) - g{x,Xm)]crmit)dt\, o -> 0. (36) 
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If k\x — Xm\ ^ 1, and fc > is fixed then 

\g{x,Xm)\ = 0{- ] , ), \g{x,t) - g{x,x^)\ = 0{- , ), 

\x-Xm\'-'^ \x-XmVl^ 

and Qm 7^ 0, so estimate ([55]) holds. 
If 

|a; — a;„i| ^ d > a, 

then ^ 

= O(t-t-), |.g(i,i) -5(i,2;m)| = O(-), 
In- a 

a 

as follows from the asymptotics of Hq^t) = O(ln-i) as r 0, and from the 

formulas = -Hl{r) = 0(i) as r ^ 0. Thus, <^ holds for |i - i™| > 

d » a. 

Consequently, the scattering problem is reduced to finding numbers Qm, 
l<m< M. 

Let us estimate Qm asymptotically, as a — >■ 0. To do this, we use the exact 
boundary condition on Sm , which yields 

-Ue{s)= / g{s,t)aj{t)dt, s€Sj. (37) 
Jsj 

The function Ue{s) is twice differentiable, so 

Ue{s) = Ue{Xj){l + 0{ka)). 

Neglecting the term 0{ka) as a — ?> 0, rewrite equation ([BT)) as 

~Ue{xj)= I g{s,t)aj{t)dt. (38) 

This equation is similar to (1191) : the role of uq{0) is played by Ue{xj). Repeating 
the argument, given in Section 2, one obtains a formula, similar to ([28]): 

_WM[i + o(l)], a->0. (39) 
m - 

a 

Formula, similar to (f^ . is 

27r 

u{x) ^ Uq{x) - T- g{x, Xrn)Ue{Xm), a 0. (40) 

m - — ' 

a m=l 

The numbers Ue{xm), 1 < w < M, in (|40p are not known. Setting x ~ Xj in 
(gOl), neglecting o(l) term, and using the definition (pij) of the effective field, 
one gets a linear algebraic system for finding numbers Ue(xm): 

Ue{Xj) ^ Uo{Xj) - j-j^ g{Xj,X„i)Ue{Xm), ^<j<M. (41) 
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This system can be solved numerically if the number M is not very large, say 
M < O(103). 

If M is very large, M = M{a) — oo, a — !■ 0, then we derive a linear integral 
equation for the limiting effective field in the medium obtained by embedding 
many cylinders. 

Passing to the limit a — > in system (|^T|) is done as in [13^. Consider a 
partition of the domain D into a union of P small squares Ap, of size b = &(a), 
6 d » a. For example, one may choose b = 0(a^/'^), d = 0(a^/^), so that 
there are many discs Dm in the square Ap. We assume that squares Ap and 
Aq do not have common interior points ii p q. Let yp be the center of Ap. 
One can also choose as yp any point Xm in a domain Dm C Ap. Since Ue is 
a continuous function, one may approximate Ue{xm) by Ue{yp), provided that 
Xm C Ap. The error of this approximation is o(l) as a — >■ 0. Let us rewrite the 
sum in (14111 as follows: 



2n 



^Jt^ 9(^J^^"^)^eiS:m) ^"^T^ X! 9(.Xj,yp)ueiyp):^ ^ 1, (42) 

a m^j p=l a x,„eAp 

and use formula ^ in the form 

^ ^ l = Ar(yp)|Ap|[l + o(l)]. (43) 

Here |Ap| is the volume of the square Ap. 
From JUl and (HSl one obtains: 



277 ^ 

9{xj,Xm)ue{xm) = 2n ^ g{xj , yp)N {yp)ue{yp)\Ap\[l + o{l)]. (44) 

a m^j p=i 

Xj ^ Ap 

The sum in the right-hand side in (j44p is the Riemannian sum for the integral 

lim y]5(*j,yp)^(yp)^^e(yp)|Ap| = / 9{x,y)N{y)u{y)dy, u{x) ^ Yun u^{x). 

(45) 

Therefore, system (|^T|) in the limit a — > yields the integral equation for the 
limiting effective field 

u(x) = uo(i) - 27r / g{x,y)N{y)u{y)dy. (46) 



D 



One obtains system (HTI) if one solves equation (|46p by a collocation method. 
Convergence of this method to the unique solution of equation (|46l) is proved in 
|10j . Existence and uniqueness of the solution to equation (|46p are proved as in 
[6], where a three-dimensional analog of this equation was studied. 
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Applying the operator A2 + to equation (H^ yields the following differ- 
ential equation for u{x): 

A2u{x) + k'^u{x) - 2ttN{x)u{x) = 0. (47) 

This is a Schrodinger-type equation, and u{x) is its scattering solution corre- 
sponding to the incident wave uq = e*'^*'. 

Let us assume that N{x) — N is a constant. One concludes from (H7|) that 
the limiting medium, obtained by embedding many perfectly conducting circular 
cylinders, has new parameter nf^ := — 2ttN. This means that k'^ = + fc| 
is replaced by k'^ := k'^ — 2ttN. The quantity fc| is not changed. One has 
P = w^n^, k'^ = uj^Uq. Consequently, ri? /tiq = (/c^ — 2nN)/k'^. Therefore, the 
new refraction coefficient is 

= n^(l - 27riV/c-2), (48) 

Since the number iV > is at our disposal, equation p8)) shows that choos- 
ing suitable N one can create a medium with a smaller, than tiq, refraction 
coefhcient. 

In practice one does not go to the limit a — ?> 0, but chooses a sufficiently 
small a. As a result, one obtains a medium with a refraction coefficient n^, 
which differs from (|48| a little, lima_>o = . 

4 Conclusions 

Asymptotic, as a — > 0, solution is given of the EM wave scattering problem by 
many perfectly conducting parallel cylinders of radius a. The equation for the 
effective field in the limiting medium obtained when a — > and the distribution 
of the embedded cylinders is given by formula ([2])- The presented theory gives 
formula (fi5)l for the refraction coefficient in the limiting medium. This formula 
shows how the distribution of the cylinders influences the refraction coefficient. 

5 Appendix 

Let us derive formulas (IT?l) - (fT5]) . Look for the solution to ©-(H]) of the form: 

Ei=e'^''-E^{x,y), E2 ^ e^'"'' E^ix.y), E^ ^ e^'''-u{x,y), (49) 

Hi=e'''''Hiix,y), H2 = e'''^' H2{x,y), i/3 = 0, (50) 
where fcs — const. Equation Q yields 

Uy - ik3E2 ^ iijJfJ-Hi, -Ux + ik^Ei = iuj^j.H2, E2,x = Ei^y, (51) 

where, e.g., Ej,x ■= Equation ^ yields 

ik^H2 ~ iujeEi, ik^Hi = —iujeE2, H2,x ~ ^i.y = —i'-^su. (52) 
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Excluding Hj, j — 1,2, from ([?T|) and using (15^ . one gets 

El = -^Ur, E2 = -^Uy, E3 = u, (53) 

iJi = H, = -^zi., iJ3 - 0. (54) 

Since Ej = E^e'-''^'' and iJj = Hje'''^'', formulas ([ll])-(IT5]) follow immediately 
from (I551)-(l54 
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